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We discuss the weak coupling expansion of a one plaquette SU{2) lattice gauge theory. We show 
that the conventional perturbative series for the partition function has a zero radius of convergence 
and is asymptotic. The average plaquette is discontinuous at = 0. However, the fact that SU{2) 
is compact provides a perturbative sum that converges toward the correct answer for positive g^. 
This alternate method amounts to introducing a specific coupling dependent field cut, that turns 
the coefficients into g-dependent quantities. Generalizing to an arbitrary field cut, we obtain a 
regular power series with a finite radius of convergence. At any order in the modified perturbative 
procedure, and for a given coupling, it is possible to find at least one (and sometimes two) values of 
the field cut that provide the exact answer. This optimal field cut can be determined approximately 
using the strong coupling expansion. This allows us to interpolate accurately between the weak 
and strong coupling regions. We discuss the extension of the method to lattice gauge theory on a 
_D-dimensional cubic lattice. 
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I. INTRODUCTION 

Lattice gauge theory incorporates essential features of 
the strong interactions at short distance (asymptotic free- 
dom) and large distance (confinement). Expansions in 
1/(3 = /2N and /3 usually provide good approxima- 
tions for the average value of gauge invariant quantities 
in the limit of small or large /3. However, calculations 
in the intermediate region often require a numerical ap- 
proach. 

There exists a general method for calculating Wilson' 
s or Polyakov's loops in powers of 1//3 0j in pure SU {N) 
gauge theories (see Ref. for a more complete set of 
references on lattice perturbation theory). 

Much effort has been devoted calculating 



P ^ ( (1/A^p) ^(1 - {llN)ReTrU^) 



(1) 



where Up denotes the usual product of links along a 1 x 1 
plaquette and Np the number of plaquettes. P can be ob- 
tained by taking the derivative with respect to /3 of the 
free energy density. Exact calculations of the coefficients 
of P up to order 3 in 1//3 and numerical calculations 
at order 8 Q and 10 's'l are available. The accuracy of 
the weak and strong coupling expansions at successive 
orders is shown in Fig. ^for S'[/(3) in 4 dimensions. The 
figure makes clear that in the region 5 < /3 < 6, none of 
the expansions (in powers of (3 or 1//3) is accurate. Un- 
fortunately, this region is precisely the "scaling window" 



*AIso at Obermann Center 
I yaniiick-meurice®uiowa.edu| 



for Advanced study; 



where one can extract information relevant to the con- 
tinuum limit. In addition, the convergence of the weak 
coupling expansion is not completely understood. The 
analysis of the numerical series may suggest the pos- 
sibility of a finite radius of convergence (the center of the 
circle being — Q). This possibility is not expected on 
general grounds and is in contradiction with the disconti- 
nuity of P when changes sign 0, [1] . We are not aware 
of any independent argument in favor of a finite radius 
of convergence, and the most likely outcome is that the 
factorial growth of the series takes over at higher order. 
In Ref. 6, this order is estimated to be approximately 25, 
which is out of reach of numerical calculations. 
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FIG. 1: P versus (3 for 5(7(3) in 4 dimensions. The solid line 
represents the numerical values; the dashed lines on the left, 
successive orders in the strong coupling expansion; the dot- 
dash lines on the right, successive orders in the weak coupling 
expansion. 
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In this article, we show that for a SU (2) lattice gauge 
model on a single plaquette, the weak coupling expan- 
sion is asymptotic (has zero radius of convergence), but 
that it is possible to modify the perturbative procedure 
in order to get a convergent sum, which is an "expan- 
sion" in powers of but with /3-dependent coefficients, 
that is accurate even in the strong coupling region. This 
work is motivated by recent results obtained in the case 
of scalar field theory ^ where the answer for similar 
questions in the case of nontrivial models can be guessed 
correctly by considering a single site integral. This is 
briefly reviewed in Sec. ^1 The main point is that the 
large order behavior of perturbation theory is related to 
large field configurations and that by cutting-off these 
configurations appropriately, we can obtain a series that 
converges to a value exponentially close to the exact one 
(for instance, errors of order e~^'^"^<^^ for the simple in- 
tegral discussed in Ref. 0]). Hereafter, we follow the 
same path for lattice gauge theories. 

The SU{2) model is introduced in Sec. IIIII where we 
also discuss its connection to Bessel functions. It is worth 
noting that for compact groups, there are no large field 
contributions and consequently, the factorial growth of 
the perturbative series comes from adding the tails of in- 
tegration, as done in asymptotic analysis of integrals 
and in the conventional procedure used in lattice per- 
turbation theory 1]. This is explained in Section Hvl In 
many quantum problems, the lack of convergence can be 
traced to the behavioiMjf the model at negative coupling 
(see, however, Ref. jl^ for a proper definition). This 
question has been discussed for lattice gauge models in 
4 dimensions 0. In Sec. El we argue that there should 
be an essential singularity at zero coupling for the one 
plaquette model. In Sec. IVII we show that the regu- 
lar perturbation series (with the integration tails added) 
misses "instanton effects" of the form (3~^e~'^^ . 

In Sec. IV 111 we propose to modify the conventional 
perturbative method by introducing a field cut. With 
this modification, the series converges toward a value 
which in general is different than the exact one. How- 
ever, at a given order in the weak coupling expansion, 
it is possible to pick an optimal field cut such that for a 
given coupling the answer is exact. For the integral stud- 
ied in this article, we found at least one solution at every 
order. This is not necessarily the case in general. For the 
integral studied in Ref. ^3|> were able to prove that 
no such solution exists at odd order and that we could 



and discuss three practical ways to calculate the modified 
coefficients. 



only minimize the error in that case. In Sec. IVIIII we 
use the strong coupling expansion to determine approxi- 
mately this optimal field cut. In this approach, the field 
cut is given as a power series in /?. A numerical study in- 
dicates that this series has a finite radius of convergence 
which increases with the order in 1//3 considered. The 
method that we propose allows us to interpolate between 
the weak and strong coupling region. This is depicted in 
Fig. 12 which is the prototype of what we expect to ac- 
complish in general. In the conclusions, we consider the 
implementation of the method for _D-dimensional models 
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FIG. 2: Number of significant digits for the SU{2) one pla- 
quette integral, at order 7 in the weak coupling (dotted line 
W7), at order 2 in the strong coupling (empty circles SC2) 
and at order 7 of the modified perturbative method proposed 
here with an optimal field cut determined pointwise using the 
strong coupling expansion at order 2 (solid line W7S2). 



II. MOTIVATIONS 

A common challenge for quantum field theorists con- 
sists in finding accurate methods in regimes where exist- 
ing expansions break down. In the renormalization group 
language, this amounts to finding acceptable interpola- 
tions for the flows in intermediate regions between flxed 
points. A discussion of this question for lattice gauge the- 
ories can be found in Refs. hj,hj. In the case of scalar 
field theory, the weak coupling expansion is unable to 
reproduce the exponential suppression of the large field 
configurations operating at strong coupling. This prob- 
lem can be cured QJ by introducing a large field cutoff 
(^max which eliminates Dyson's instability. One is then 
considering a slightly different problem, however a judi- 
cious choice of 4>max can be used to reduce or eliminate 
^3 the discrepancy with the original problem (i. e., the 
problem with no field cutoff). This optimization proce- 
dure can be approximately performed using the strong 
coupling expansion and naturally bridges the gap be- 
tween the weak and strong coupling expansions. 

The study of the simple integral 
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provides a good understanding about the role of large 
field configurations in the perturbative series. It helps 
identifying general features of the the effect of a field 



3 



cut. In particular, the dependence of the accuracy of 
the modified perturbative series on the couphng and the 
field cut, is qualitatively very similar for the integral, 
the anharmonic oscillator and the hierarchical model in 
3 dimensions (see the similarity among the three parts of 
Fig. 2 in 9). 

In order to generalize this procedure to gauge theory, 
we will first consider the simplest possible gauge model, 
namely the one plaquette integral 

Z{p,N) = / J]d[/ie-'3(i-^«-r^^^) , (3) 

After fixing the gauge so that U = 1 on three sides of the 
plaquette, Z becomes an integral over a single link 

Z(/3, N) = j df/e-/3(i-^fl-i^-^) (4) 

For an arbitrary gauge fixing prescription, TrU becomes 
Tr{U{g)U) with g arbitrary and Z is (^-independent by 
virtue of the invariance of the Haar measure dU. This 
integral and its moments appear in the strong coupling 
expansion _^13, 15, 16, 17] and in the mean field treat- 
ment [13 ofsU{N) gauge theories. In the one plaquette 
model, 

P = --il..Z. (5) 

The accuracy of successive orders in the f3 and 1//3 de- 
scribed in the following sections is shown in Fig. |3|that 
can be compared with Fig. ^ 



III. THE MODEL CONSIDERED HERE 

In the following, we specialize the discussion to the 
case N — 2 for which the Haar measure is very simple. 
From now on, the reference to iV will be dropped and we 
will use the notation Z{f3) for Z{(3, 2). The explicit form 
is: 

1 r'^'^ 

Z{I3) = - dujsm\uj/2)e-^'^^-^°'^^/^^K (6) 
I" Jo 

Setting u = cos(aj/2), 

Z{f3) = - C du^\ - u2e-'5(i-") (7) 

J-1 

and one recognizes from Eq. (8.431) of Ref. 0| that the 
integral can be expressed in terms of the modified Bessel 
function Ii. 

Z(/3) = 2c-^/i(/3)//3 (8) 

Using the Taylor expansion Eq. (8.445) in Ref. we 
can write 

As in the case of the integral of Eq. (0) , the presence of 
the factorial at the denominator implies that the strong 
coupling expansion (in powers of P — 4/3^) converges 
over the entire complex plane. 

IV. THE WEAK COUPLING EXPANSION 
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FIG. 3: P versus /3 for SU{2) on one plaquette. The solid line 
represents the numerical values; the dashed lines on the left, 
successive orders in the strong couphng expansion; the dot- 
dash lines on the right, successive order in the weak coupling 
expansion. 



Assuming /3 > 0, we set t = /?(1 — it) in Eq. Q yields 
Z(/3) = (2//3)3/2- / dtt'/^e-'^l-{t/2P) (10) 



If we expand the square root in the integral and exchange 
the sum and the integral (the validity of this procedure 
will be discusssed in Sec. IVII|) . we obtain a converging 
sum: 

Z{(3) = (/3^)-3/22i/2 ^ M2P)P-' , (11) 



1=0 



with 



The convergence of the sum in Eq. (|ll|l can be estab- 
lished from the bounds 
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Z + 3/2^ ^ ^ io Z-f-3/2 



(2/3) 



;+3/2 



(13) 
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and the fact that r(; + 1/2)//! < 1 for / > 1. (This is 
true for I = 1 and can be proved by induction muhiplying 
the inequahty by {I + l/2)/(Z + 1) < 1). Consequently, 
the sum converges at the same rate as ^ Z^^. Note that 
as in the case of the ground state of the quantum me- 
chanical double well, the first term is positive but all the 
remaining terms are negative. 

Obviously, Eq. (|ll|l is not a power series in since 
the "coefficients" A; (2/3) are /3-dependent. We can now 
obtain the conventional asymptotic expansion by two dif- 
ferent methods. The first one consists in adding the tails 
to the integrals in Eq. (|12|) , or in other words by replacing 
the incomplete gamma function by the gamma function. 
This is a standard procedure in asymptotic expansions of 
integrals pH . 

On then obtain the asymptotic expansion 

Z{P) ~ (/3^)-3/22i/2 X 

The terms of this sum now grow like Z!/2' and the series 
is asymptotic. As all the signs are negative for / > 1, 
the Borel transform has singularities on the positive real 
axis. 

It is instructive to rederive the expansion of Eq. H14(l 
by following the steps of lattice perturbation theory [l| . 
We first set oj — gA in Eq. ^ and expand the action and 
the Haar measure in powers of g. This leaves us with the 
integral of a power series in g over the range to 2T:/g 
for A. The asymptotic series ()14f) is then recovered by 
letting the range of integration go to infinity. As the two 
methods amount to calculate the coefficients with differ- 
ent variables of integration, we obtain the same series, as 
can be checked explicitly up to high order. We emphasize 
that in lattice gauge theory with compact groups, there 
are no large field contributions. It is only for practical 
reasons that the tails of integration are added. In the one 
plaquette example, calculating A; (2/3) instead of Ai{oo) 
is a very minor problem, however, this is a technical chal- 
lenge in the case on a D-dimensional lattice 

V. BEHAVIOR AT NEGATIVE f3 

From the integral representation Eq. the change 
/3 — *■ —(3 can be made by changing u — > — u and multi- 
plying by e^^. This implies, 

Z{-f3) - e'^'ZiP) , (15) 

and 

P(/3) + P(-/3) = 2 (16) 

A similar equation can be found for a SU{2) pure 
gauge model on a cubic lattice. Since ZiTO/3^+oo-P(/3) = 0, 
the limits g"^ — > 0^ differ by 2 and a converging series in 
g about is impossible. 



The discontinuity in the values of P near g^ = 
appears in a simpler model where the integration over 
SU{2) is replaced by a sum over the two elements of its 
center: 

= J2 e-^(i-i«-^^'-^) = 1 + e-'^ . (17) 
u=±t 

This implies 

2 

Pcenter = Y+'^ ' '^^^'^ 

The center model satisfies Eqs. (|15|l and (|16() . Note 
that Eq. I|17|l makes clear that Zcenter has an essential 
singularity at 5 = 0. The asymptotic behavior of Pcenter 
at large |/3| in the complex plane is 2(1 — e^^ -I- . . . ) if 
Ref3 < 0, and 2e"^^ -I- . . . , if i?e/3 > 0, with Stokes lines 
running along the imaginary axis. 

This simplified example makes clear that the usual per- 
turbation series is obtained by making modifications of 
order e"^/^ (the effect of the tails of integration). We 
now proceed to estimate the order e~^^ corrections to 
the integral over the whole SU{2). 

VI. ACCURACY OF REGULAR 
PERTURBATION THEORY 

In the study of scalar models , we have shown that if 
we plot the accuracy of perturbative series at successive 
orders, an envelope setting the boundary of the accu- 
racy that can be reached at any order using the usual 
perturbation theory appears. In the case of the quan- 
tum mechanical double-well, this envelope coincides very 
precisely with the instanton effect. We expect the lim- 
itation in accuracy to be of the general form g^e~^/^ . 
For the model considered here, the limitation of accuracy 
has this generic form and we will see that the effect is of 
order (3~^e~^^ . 

For (3 not too small, the low orders the asymptotic 
series Eq. (|14f) overestimate Z. As we let the order in- 
crease, the series crosses the true value and then start 
to grossly underestimate the true value. At each order, 
there is a special value of (3 for which the truncated se- 
ries coincides with the exact answer. This explains the 
"spikes" seen in Fig. ^ 

If we assume that for a particular value of /3, the con- 
verging sum, Eq. Hll|) with the integrals running from 
to 2/3, truncated at order K is a. good approximation of 
Z{P), then the error 6Z{P, K) made by using the regular 
perturbative series, Eq. H14() with the integrals running 
from to CX3, truncated at the same order, is in good 
approximation 

(5Z(/3,if) ~ (/37r)-=^/22i/2 ^ 
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FIG. 4: Number of correct significant digits as a function 
of /3 at successive orders of the regular perturbative series 
Eq. (inj for Z{f3). As the order increases from 1 to 15, 
the curves {Wl, W2, . . . ) get Ughter. The thick sohd hue is 

Integrating by parts, dropping terms of order /?~^ and 
summing the resulting series, we obtain 

SZ{I3, K) « ^^6-2^/3-12^-3/2 , (20) 

with 

1=0 

The coefficient Ak slowly decreases when K increases. 
For instance, A5 = 0.872 . . . , Aw = 0.624 .... In the 
limit of large K, Ak becomes zero. This comes from the 
resummation 

V = r dtt-^'\l - c-*) = 2^1/2 , (22) 

which is exactly the I = term. In practice, when the 
order is not too large, fi~^e~'^^ is a good order of mag- 
nitude estimate for the envelope discussed above as can 
be seen in Fig. 0] 



VII. MODIFIED MODEL WITH AN 
ARBITRARY FIELD CUT 

In this section,we consider a modified partition func- 
tion Z{P,t^ax) where the integration range of Eq. 1(1171) 
takes a fixed, /3-independent, value tmax- When tmax < 
2/3, the Taylor series for the square root converges abso- 
lutely and uniformly over the whole range of integration. 
It is thus justified to interchange the sum and the integral 



and we have 

00 

Z{f3,t^ax) = {p7r)-^^^2^^^J2Mtmax)r' ■ (23) 

1=0 

The original partition function as expressed in Eq. Hll() 
is obtained by taking the limit tj^ax ~* 2/3. Since the in- 
tegral with upper boundary tmax is obviously continuous 
in that limit and since we can use the suppression 
to prove the continuity of the sum, the validity of Eq. 
(|^ extends which justifies Eq. (fTT|l . If t^ax > 2/3, 
the sum diverges and the integral is ill-defined. The reg- 
ular perturbation series is obtained by taking the limit 
tmax oo. In the graphs, we use the notation "WK" for 
the truncation of the regular perturbative series at order 
K. One should however keep in mind that, for instance 
in W7, the last term is of order (l//3)''+3/2 due to the 
prefactor /3^3/2 Eq. if^ . 

Eq. H23|) is now a regular series in (1//3). It has a 
finite radius of convergence. In order to calculate this 
radius, we notice that for large /, /q™°^ dtt' gets most of 
its contribution from the region between tmax{^ — 1/0 
and tmax- Consequently, one can replace e"*ti/2 j^y 

1/2 

e~*"""'tmax without affecting the asymptotic behavior 
of the coefficients of the series. If we perform this 
change directly in the integral Eq. H10|l . the integral 
can be calculated explicitly. One can then conclude 
that Z [(3^ tmax) has a non-analytical part proportional 
to (1 - (Uax/2/3))3/2. The series defined by Eq. ^ 
converges if (1//3) < [2 /tmax)- Numerical studies of 
the series with conventional estimators confirm this ar- 
gument. Note that the finite radius of convergence of the 
series Eq. (refeq:mod) is not in disagreement with the 
discontinuity of the original model at 1//3 = 0, because 
this series coincides with the original model only when 

(1//3) = {2/tmax)- 

Can a truncation of the series of Eq. (|23|) at order K 
be a good approximation of the original integral Eq. ^1 
The answer depends on iC, tmax and (3. It is clear that if 
K is large enough and (3 shghtly above 2/ tmax ^ then one 
should get a reasonable approximation. This statement 
is confirmed by Fig. |31 where the accuracy of Eq. I|23|) 
with tmax — 8 truncated at orders 1 to 15 is displayed 
as a function of /3. In this particular case, the values of 
/3 > 4 are within the radius of convergence. As the order 
increases, the spikes in this region (the right half of the 
Fig. [SJ move toward 4. In addition, there is another set 
of spikes, outside the radius of convergence (on the left 
half of the Fig. |SJ| and moving in the opposite direction 
when the order increases. 

A more global information regarding the location of the 
spikes is displayed in Fig. |H| It shows that the "second 
solution", outside the radius of convergence, disappears 
beyond some critical value of (3. As the order in the 
weak coupling increases, both solutions get closer to the 
tmax = 2/3 line. 
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FIG. 5: Number of correct significant digits as a function of 
/3 for a fixed cut tmax ~ 8. As the order increases from 1 to 
15 {Wl, W2, . . . ), the curves become fighter. 




FIG. 6: Location of the exact matching between the series 
Eq. 11231 at order 6, 7 and 8 and Z{P) in the f3-tmax plane. 
The dashed lines represent the solution within the radius of 
convergence and the empty circles the other solution. 

VIII. OPTIMIZATION 

In this section, we discuss an approximate method to 
find the optimal value of t^ax corresponding to a given 
order K and a given value of f3. In a general situation, 
we do not know accurately the value of the quantity that 
we are calculating (the equivalent of Z here). Conse- 
quently, we need to find an approximation that allows 
us to consistently estimate this quantity and the way its 
order K approximation changes with the field cutoff in 
order to impose an approximate matching condition. For 
this purpose, we will use the strong coupling expansion 
(power series in /3) which provides information comple- 
mentary to the weak coupling. Now, the crucial point is 



that the field cut allows us to control the (1//3) in the 
integral Eq. (|10|l , because (except in the exponential) all 
the factor {1//3) appear together with a factor t. In other 
words, except for the exponential, it is a function of t/ (3. 

We would like to match the the strong coupling expan- 
sion 

Z{(3)^l-P+{5/8)(3' + ... (24) 

discussed in Sec. IIIII with the truncated expansion of 
Eq. (ISSJ which can be rewritten as 

The control oi s — t/ P can be achieved by imposing that 
imax/ P is approximately constant. We can then improve 
order by order in [3 by setting 

{t^ax/P) = c^{K) + c^{K)(3 + ... (25) 

The only non-trivial part is to solve the zeroth order (in 
P) equation 

Fk{co{K)) = 1 , (26) 

with 

We have checked that for K going from 1 to 40, Eq. 1)26(1 
has exactly two solutions on the positive real axis with 
one solution on each side of 2. As K increases, the 2 roots 
get closer. They should coalesce at 2 in the large K limit. 
This follows from the fact that Foo(2) = 1 and F'^{2) = 
(as can be shown by using the the same method as for Eq. 
((22|l V The higher order coefficients ci{K) corresponding 
to each of the two solutions at order can then be found 
by solving linear equations. This procedure provides an 
approximate value of the optimal tmax which apparently 
converges toward the correct numerical value. This is 
illustrated in Fig. [T] If we plug the two approximate 
values of tmax of Eq. H25|l in Eq. H23|l truncated at order 
we obtain approximate values of Z{P). The accuracy 
of this procedure is displayed in Fig. |Slin the case K = Q. 
It appears clearly that the first solution (the one within 
the radius of convergence with tmax < 2/3) is significantly 
more accurate than the other solution (with tmax > 2/3). 
Similar features were observed for K up to 20. 

We can now compare the accuracy of the method pro- 
posed here with the weak and strong coupling expan- 
sions. The case K — Q is, displayed in Fig. |(5J. In the 
weak coupling region (/? > 10) the accuracy of our pro- 
cedure merges with the regular perturbation series. In 
the strong coupling region (/? < 0.1), our procedure is 
more accurate than the regular expansion in powers of P 
by several significant digits. As /3 — > 0, the accuracy of 
our procedure with tmax/P determined at order m in P 
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FIG. 7: Approximate locations in the (/3, imai) plane of the 
matching between the order 6 weak coupling expansion and 
Z{/3). The two solid lines are the two numerical solutions 
at that order (as in Fig. The dash line (empty circles) 
represent the first (second) approximate solutions at order 
0, . . . , 4 in /3. 



FIG. 9: Significant digits obtained from Eq. 1231 truncated 
at order 6 using the first solution for tmax/P at order to 3 
compared to the weak coupling expansion at order 6 (dotted 
line W6) and the strong coupling expansion at order to 2 
(empty circles SC) 



IX. ASYMPTOTIC BEHAVIOR OF ci{K) 
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FIG. 8: Significant digits obtained from Eq. 12311 truncated 
at order 6 using tmax/P at order (squares), I (circles), 2 
(triangles) and 3 (stars). The first solution with tmax < 2/3 
is showed with filled symbols, while the second solution is 
showed with empty symbols. 



In this section, we study empirically the asymptotic 
behavior of the coefficients ci{K) appearing in the ex- 
pansion of tmax/P Eq. (123) • At fixed K large I , Fig. ^ 
suggests that 



ciiK) cx iG{K)y 



(28) 



In addition, it appears that G{K) decreases with K ap- 
proximately like 1/K. This behavior implies a finite 
radius convergence G{K)^^ for the P expansion in Eq. 
(|25|1 . increasing linearly with K. This is good news for 
the interpolation between the weak and strong coupling 
region since as we increase the weak order K, we increase 
the range of validity in /?. 

The large- if behavior of ci{K) has also been studied 
numerically. The results for Z up to 5 are shown in Fig. 
rm At fixed I large K , the linear fits used in Fig. 
suggest that 



co{K)^2 + 0{l/K) 



(29) 



and 



increases at the same rate as the regular strong coupling 
expansion at order m in /3 maintaining the difference in 
accuracy approximately constant. In the intermediate 
region where none of the conventional expansions work 
well (except at the perturbative spike), our procedure 
maintains a very good accuracy interpolating smoothly 
between the two regimes. 



Q(if) ocX-'-i+"(') , (30) 

for I > 0, with a{l) small. This behavior is expected, 
since as the order increases, we are getting close to the 
exact expansion Eq. Hll|) with tmax = 2/3 (co — 2, ci — 
for Z > 0). The values oia{l) decrease when we reduce the 
set of points fitted to larger values of K. If we use K — 35 
to 45 for the fit, we have approximately a{l) ~ I/IO. 
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FIG. 10: In \ci{K)\ versus / for if = 2, 6 , 10. 
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FIG. 11: In \co(K) - 2| versus \nK; in |ci(/C)| versus InTf for 
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X. CONCLUSIONS 

We have shown that for the one-plaquette model, the 
introduction of a properly chosen field cut can provide 
a high accuracy in regions where the usual perturbative 
method is not accurate. The strong coupling expansion 
provides an efficient way to determine the optimal cut 
and interpolate between the small and the large (3 re- 
gions. Apparently, the accuracy of the calculations im- 
proves whenever we increase the order in either the weak 
or the strong coupling expansion. Given these positive 
results, we are compelled to implement the method in 
the case of lattice gauge theory on a D dimensional cu- 
bic lattice. Two steps are necessary. First, we need to 
define the theory with a field cut (the analog of Eq. 1(1^ 
with 2/3 replaced by tmax)- Second, we need to expand 
relevant quantities such as P for the modified theory in 



powers of 1//3 (the analog of Eq. H23|) '). Note that in the 
calculation of P using a perturbative series, the complex 
zeroes of Z will play an important role. This question 
remains to be examined in detail. 

The implementation of the first step is straightforward. 
One can insert 1 in the partition function in the following 
way: 



1-n 



dtp5{l - {l/N)ReTr(Up 



(31) 



If we could perform the integration over the Uunk, we 
would get an effective theory for the new variables tp. 
Note that the procedure is gauge invariant since TrlJp is. 
The "size" of a configuration can be defined in several 
ways. For instance, we could use the value of MaXp {tp} 
or {1/Np) J2p^p ^'^ decide if we have a large or a small 
field configuration. We can then order the configurations 
according to the chosen indicator. Given a (sufRciently 
large) set of Monte Carlo configurations, one can define 
the expectation values with a cut by averaging only over 
configurations for which the chosen indicator is below 
a certain value. The correlations between the two size 
indicators mentioned above are now being studied for 
SU{3) in 4 dimensions. 

The implementation of the second step requires tech- 
nologies that are now being developed in the scalar case. 
As it seems only possible to make analytical calculations 
for small or large field cuts, numerical methods seem un- 
avoidable. For the purpose of independent verification, 
it is important to consider different methods. We are 
presently working on three different approaches: 

1 . The conventional approach but with the A"^ hav- 
ing a finite range of integration. This type of ap- 
proach works well in the scalar case 20j 

2. The stochastic approach 4] where A'^ is expanded 
as power series in 1/(3. For the lowest order field, 
the implementation of a cut is obvious but not for 
higher order fields. This problem is being consid- 
ered with simple examples. 

3. Fits from numerical data at large (3. This method 
[2ll | allowed to extract at least 2 coefficients of con- 
ventional perturbation. As we mentioned above, 
it is easy to implement the field cut with Monte 
Carlo methods. The advantage of this method is 
that it does not require the use of the Campbell- 
Baker- Hausdorff (for a short review see Ref. [221] ') 
formula. 

We expect that the use of theses three methods will al- 
low us to contruct perturbative series with a finite radius 
of convergence as above. We hope that this radius of 
convergence will be sufficiently large to reach the scaling 
window. Ultimately, we expect to be able to replace the 
numerical calculation of the coefficients by approximate 
analytical formulas, as it seems possible to do in the case 
of the anharnionic oscillator |2^ . 
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